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Abstract 

We study the fluctuations of a random surface in a stochastic 
growth model on a system of interlacing particles placed on a two di- 
mensional lattice. There are two different types of particles, one with 
a low jump rate and the other with a high jump rate. In the large 
time limit, the random surface has a deterministic shape. Due to the 
different jump rates, the limit shape and the domain on which it is 
defined are not smooth. The main result is that the fluctuations of the 
random surface are governed by the Gaussian free field. 

1 Introduction 

The Gaussian free field is commonly assumed to be a universal field describ- 
ing the fluctuations of random surfaces appearing in a wide class of models 
in statistical physics. However, rigorous proofs are only known for some par- 
ticular integrable models. For example, interesting progress has been made 
on dimer models on bipartite planar graphs (see [15] for a survey and a list 
of references). 

In the present paper we study a random surface appearing in a stochastic 
growth model on a system of interlacing particles and show that its fluctua- 
tions are governed by the Gaussian free field. This work is partially inspired 
by [8], where the authors introduced a general model in 2+1 dimensions that 
connects the random surfaces of the type that occur in the dimer models, 
with the random growth of particle systems on a one dimensional lattice 
(e.g. exclusion processes). Here we will discuss a particular specialization 
of that model. 
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Figure 1: The left figure shows the initial configuration of the particles. The 
right figure shows a possible configuration after a short time. The dashed 
line separates the particles with a higher jump rate (above) from the particles 
with a lower jump rate (below). 

In the initial configuration the particles are placed on the grid Z x N 
and are densely packed in a triangular region as shown in the left picture 
in Figure [TJ Each particle has an exponential clock and when the clock 
rings, the particle attempts to jump to the right by 1. To ensure that 
the interlacing of the particles on subsequent horizontal levels is preserved, 
each attempted jump is subject to certain rules (specified later on, but 
in short: each particle is blocked by the particles below, but pushes the 
particles above). The central feature of our model is the fact that we have 
two different types of particles: slow and fast. More precisely, we draw a 
horizontal separating line and equip the particles above that line with a 
higher jump rate than the particles below. As a consequence the particles 
above the separating line tend to move faster. See [8j for a thorough analysis 
for the situation in which all particles have the same jump rate. 

For large time, the particles will be distributed on a large domain that 
is shown in Figure |4} See also Figures [2] and [3] for sample configurations. 
There exists a critical time, after which the limiting domain develops a cusp. 
This means that after the critical time, there is a group of fast particles that 
is drifting away from the slower particles below the separating line, while 
other fast particles are held back by the slow particles below due to the jump 
rules. The latter situation is illustrated in Figure [3] and the right picture in 
Figure |4| 

Our main results on the long time behavior are expressed in terms of 
a height function that integrates the particle configuration. The graph of 
the height function defines the random surface that has our interest. We 
show that this height function has a deterministic shape in the large time 
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limit. The difference in speed for the two types of particles induces a jump 
discontinuity in the normal to the limit shape at the line separating the slow 
and fast particles. 

In order to describe the limit shape, we introduce a bijection that maps 
the limiting domain to the upper half of the complex plane. This map 
constitutes a natural complex structure on the system. Away from the 
separating line it satisfies the complex Burgers equation (for more details 
on the connection between this equation and limit shapes see [16J). Due to 
the different jump rates, the bijection in our case is only homeomorphism but 
not a diffeomorphism (in contrast to the situations in for example [U I14j). 
The main result of the paper is that the fluctuation of the pushforward of 
the random surface under this map are governed by the Gaussian free field 
on the upper half plane. 

The limit shape is subject to a PDE. Away from the separating line, 
this PDE brings our model locally in the 2D anisotropic KPZ class [2] (see 
also [8] and the references therein). This constitutes an important class of 
stochastic growth models. In [28J, the author used non-rigorous arguments 
to show that the fluctuations in these models are similar to the fluctuations 
of random surfaces in the Edwards- Wilkinson class, which predicts the log- 
arithmic behavior of the height fluctuations. However, it does not predict 
the full result on the Gaussian Free field, especially in relation with the 
complex structure. The only model of this type (that the author is aware 
of) for which the Gaussian free field in the fluctuations is rigorously proved, 
is the one jump rate situation analyzed in [8]. In our situation, the PDE 
has a jump discontinuity on the separating line. It is not a priori clear if 
and how this effects the fluctuations. For example, it is not obvious at first 
(and perhaps slightly surprising), that the correlations of the fluctuations of 
the height function at a point in the slower part and at a point in the faster 
part that starts drifting away, are essentially of the same type as for points 
that are both in the slow part. The essence of our main result is that in 
all cases these correlations are given by the pullback of the Green's function 
for the Laplace operator on the upper half plane with Dirichlet boundary 
conditions by the complex structure on the system. 

Finally, we remark that in for example [BJ [T3"l PH} [15] , the presence of 
the Gaussian free field is established by computing the limiting behavior for 
the moments of the height fluctuations at multiple points. In this paper, we 
follow an alternative approach that exploits the determinantal structure of 
the process. 
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Figure 2: A random particle configuration before the critical time. The 
difference in jump rates is exaggerated for illustration purposes. 

Acknowledgements I thank Alexei Borodin for drawing my attention to 
the subject of this paper and for many fruitful discussions. 

2 Statement of results 

In this section we will state our main results. The proofs can be found in 
Sections [HE 

2.1 The model 

Let us start with describing the model. We consider dynamics on a system 
of interlacing particles placed on the grid Zxff. At each point in time, 
there are m particles on the horizontal section (-,m) for m £ N. We denote 
the horizontal coordinate of the fe-th particle (counting from left to right) 
on the m-th level by x™ for k = 1, .. .,m. In the initial configuration at 
time t = the positions are 

(2.1) x%(Q) = k-m-l, for k = l,...m, 

as shown in the left picture of Figure [TJ Each particle has an exponential 
clock and when its clock rings, it attempts to jump to the right by one. 
However, we enforce the following interlacing condition to hold at all times 

(2.2) zjffi 1 . 

To ensure that this interlacing condition holds, we impose the following 
rules. If the exponential clock of the particles positioned at rings then 
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1. it remains put if xT? 1 = xV? + 1. 

2. it jumps to the right by one and so do all particles with horizontal 
coordinate x^ 1 = x™ for I = 1, 2, . . .. 

Hence a particle is blocked by particles that are below, but it pushes particles 
that lie above. 

It remains to set the rate for the exponential clocks. The first mo hori- 
zontal levels have an exponential clock with rate 1, whereas the particles at 
the higher levels have rate 2. Hence the particles at xV? with m > tuq will 
attempt to move faster in time. 

We are interested in the long time behavior of the system. 

Remark 2.1. The above model can also be established as a tiling of the 
half plane by lozenges [5]. 

2.2 Long time behavior 

We proceed with an informal discussion on the long time behavior. Precise 
statements are formulated in the next paragraph. 

In Figures [2] and [3] we give two particular sample configurations. As these 
figures suggest, for large time the particles fill a domain as given in Figure |4j 
The slanted dashed line in Figure [4] represents the line x + m = 0. Note 
that at all times, the particles are positioned at the right of that line. The 
horizontal dashed line is the line m = mo that separates the slow and fast 
particles. As a consequence to the different speeds of the particles, there is a 
transition in the shape of the domain and the density as one crosses the line 
m = mo- The particles above that line tend to move faster than the particles 
below. Due to the interlacing condition, the particles more to the left in the 
upper part of the system are held back by the particles in the lower part. 
The particles to the right move independently from the lower part. There 
is a critical time, after which particles more to the right in the upper part 
start forming a group that starts drifting away and a cusp appears in the 
domain. This is illustrated in the domain at the right in Figure |4j 

Although the focus to this paper is on the global scale and its fluctua- 
tions, it is also possible retrieve interesting universality classes at the local 
scale. As the proofs are standard and the results are not relevant to this 
paper, we content ourselves with a brief discussion. For points in the bulk, 
the local correlation are described by extensions of the discrete sine kernel 
that fall into the class described in [3]. Near the edge, we obtain extensions 
of the Airy process (see [22] and [9] for a review). The local correlations 
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Figure 3: A random particle configuration after the critical time. The dif- 
ference in jump rates is exaggerated for illustration purposes. 



near the cusp are determined by the Pearcey process [U El El [211 12Z] • At 
the points where the domain touches the line x + m = 0, the horizontal axis 
and the line m = mo (after the critical time), the local process is governed 
by the GUE minor process [711121121!]. 

2.3 The limit shape 

To each particle configuration we assign a surface, which is the graph of the 
height function given in the following definition. 

Definition 2.2. We define the height function h : Z x N — > Nq by 



In other words, h(x,m) counts the number of particles that are at the right 
of (x, m) (including the possible particle at (x, m) ). 

For each configuration, the height function defines a stepped surface. 
Our first result is an explicit description of the limit shape of the mean E/i. 
To this end, we need the following definitions. 

Definition 2.3. Denote the upper half plane by H + = {z G C | lmz > 0}. 

Define the function F : M + — > C by 



tz + fi logO - 1) + O - hq) \og(z - 2) - (£ + //) log z, n > no, 



(2.3) 



h(x,m) = #{xf | x^ > x). 



(2.4) F(*|£,/x) 



M < Mo 





Figure 4: Two possible limiting domains T>, before (1) and after (r) the 
critical time. In the right figure, the particles in the upper part that are not 
blocked by the particles in the lower part are drifting away from the lower 
part. 



for z G H_|_. To every (£,/x) there exists at most one O G HL|_ such that 

F'(n | /z,o = o. 

If it exists we denote it by 0(£,/x) and we define 

(2.5) V = {(£, /x) | n(f , /x) G H + existe}. 

The fact that to every (£, //) there exists at most one f2 G H+ such that 
| /U,£) = 0, follows immediately after observing that F'(U \ £, /x) = is 
equivalent to a quadratic (in case zx < /xq) or cubic (in case fi> fio) equation 
in O with real coefficients. Moreover, we have the following. 

Proposition 2.4. The map (£, /x) h-> 0(£, /x) is a homeomorphism from T> 
to H + . 

Remark 2.5. As in [8l eq (1.13)] we note that that map Q, : (£, /x) h-» /x) 
solves the complex Burgers equation [16] 

{ ' (a - ft) dfi d£ ' 

where a = 1 if fi < /xo and a = 2 is /x > /xo- The jump discontinuity in a is 
due to the different jump rates. 

Remark 2.6. The boundary dV is the set of all (£,/x) sucn that ^(£>A*) is 
a double critical point of -F. For real values of £ and /x, the only possible 
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double critical points are real. This gives a way of explicitly computing the 
boundary. For SI G M, we solve the system F'(Q) = F"(Q) = for (£, /i). 
This method was used to draw the pictures in Figure [4j Moreover we have 
the following identification of points. The points Q = and 1 correspond 
to the points dT> touching the line x + m = and the horizontal axis. If 
we are beyond the critical time, we also have that $7 = 2 corresponds to the 
boundary point touching the line m = too in the cloud starts drifting away. 
The cusp corresponds to a point f2 G (1,2). At the critical time (the birth 
of the cusp), the cusp and the lowest point in the cloud drifting away meet 
and correspond to f2 = 2. 

In the following theorem we present our main result on the limit shape. 

Theorem 2.7. Let h be as defined in (|2.3|) and set 



(2.7) 



t = [Lr] 
m = [Lfjb] 
to = [Lfj, ] 
x = [L$ 



Then we have that 



(2.8) 

for (^)eD. 



lim —Kh(x,m) 



7T 



ImF((l(e,/i))=: hfon), 



The limiting height function is a continuous function of the scaled vari- 
ables. However, the normal to the surface constructed out of its graphs is 
not, as can be seen from the following result. 

Proposition 2.8. Let n = (ni, 712,713) be the norm al to the limiting surface 
defined by the graph of h defined in Theorem 2.1 Normalize n such that 
713 = 1. Then we have 



(2.J 



where 9j are the angles of the triangle formed by 

1. 0, 1 and /j) if fj, < no, 

2. 0, 2 and /j) if \x > fJ>o- 







( \ 




l-l 




w 
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Figure 5: The angles 8j in Proposition 2.8 



See also Figure [5| 

The discontinuity in the limit shape can also be seen from the PDE that 
it satisfies. See also Remark 12.51 



Proposition 2.9. The limiting mean height function h from Theorem 2.1 
satisfies the PDE 



(2.10) 

where 
(2.11) 



dh 



f(x,y) 



(f(Vh), h<hq 

sin7ra;sin7r(y — x 
it sin ny 



One easily computes that the signature of the Hessian of / is (1,-1). 
This implies that our model is locally in the class of growth models described 
by the anisotropic KPZ equation [2] as mentioned in the Introduction. How- 
ever, the PDE has a jump discontinuity at the separating line. 



2.4 Gaussian free field 

The main result of the paper is on the fluctuation of h around E/i, which 
turn out to be governed by the Gaussian free field. For a survey on the 
Gaussian free field see |23j . 

Let us first recall some definitions. Denote the Laplace operator on H+ 
with Dirichlet boundary conditions by A and consider the Sobolev space 
Wo on H + defined as the completion of the space of smooth functions with 
compact support in H + equipped with inner product 



(2.12) 



V01 • Vc/>2- 



H+ 
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Note that by integrating by parts and using the fact we have Dirichlet bound- 
ary conditions we have 



E 4 



for sufficiently smooth (p2- 

The Gaussian free field on H_|_ is a collection of centered Gaussian random 
variables {(F, 0}v}<^eW o (indexed by the Sobolev space Wo) and covariance 



(2.13) 



E[{F,<l> 1 )v{F,(fo)v] = ((h.,<f>. 



2 jV- 



The main result of the paper is that the fluctuations h — E/t in the large 
L limit are governed by the Gaussian free field. More precisely, the push 
forward of h — Kh under the map (£, fi) t-t fi) converges to the Gaussian 
free field. The proof that we present here, differs from the usual approach to 
compute the moments of the height fluctuations at different points. Instead, 



inspired by (2.13) we define the pairing of the height function with a test 



function and compute the limiting behavior of the characteristic function. 
Care should be taken here, since the function h — E/t is defined on a discrete 
set. We therefore restrict ourselves to test functions <f> that are C 2 and have 
compact support in H+. Then we define the pairing between h — E/t and 



4> by a discretization of the Sobolev inner product (2.12) and including the 
pushforward with 

For a function F:ZxN->M and a C 2 function cp with compact support 
+ , define the pairing (F, 4>) by 



m 



(2.14) (F, i 



'7T 



V F(x, m)A4> (fl(x/L,m/L)) J(x/L,m/L), 

(x,m)£LV 

where J stands for the Jacobian of the map (£,/x) h-» fi(£,/x), i.e. 

/ dRcQ dReQ s 

det 



(2.15) 



dlmQ dlmQ 



The main result of this paper is formulated in the following theorem. 
Theorem 2.10. Let cj) be a C 2 function with compact support on M + and h 



be the height function defined in (2.3). Then with (•, •) as in (2.14) we have 



(2.16) lim E[expit(h-Eh,</))] = e~^^ . 

L— >oo 

Hence, as L — )• oo, the pushfoward of the random surface defined by h — E/i 
by the map Q : T> — > H+ converges to the Gaussian free field on H+ . 
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The key fact that we use to prove Theorem 2.10 is that the process on 
the interlacing particles at time t defines a determinantal point process, as 
we will discuss in the next paragraph. 



Remark 2.11. The fact that (2.16) implies (2.13) follows by the polariza- 



tion identity for norms in Hilbert spaces and the fact that the pairing in 



(2.14) is linear. See also [23, Prop. 2.13]. 



Remark 2.12. Note that in Theorem |2.7| the height function grows with 
L as L — > oo, whereas the fluctuations in Theorem 2.10 are not scaled 
with L at all. However, the Gaussian free field is a probability measure on 
generalized functions. The pointwise limit in the scaled variables of h — Eh 
is ill-defined. Indeed, the variance at a given point grows logarithmically 
with L. In contrast, the correlation between two separate points remains 
bounded and is expressed in terms of the Green's function for the Laplace 
operator on H+ with Dirichlet boundary conditions. See also |23j. 

2.5 Determinantal point processes 

For a discrete set X a determinantal point process on X is a probability 
measure on 2 X for which there exists a kernel K : X X X — > C such that 

(2.17) Prob{X G 2 X \ Y C X) = det K(x, y) x , y£ Y 

for all finite sets Y C X. For more details on determinantal point processes 
we refer to [H Ell IH1 113 ESI EH ES]. A determinantal point process is 
completely determined by its kernel. 

The fact of the matter is that the growth model on the interlacing par- 
ticles system at time t, defines a determinantal point process on Z x N with 
kernel given by (see [8]) 



(2.18) K(xi,mi,X2,m 2 ) 

Here p m is defined by 

(2.19) p m (z) = 



Xm 1 <m 2 I P mi {w) dw 



+ 



2vri J ro p m2 (w) wxi+mi-xz-ma+l 
/ / e tw p mi (w)z X2+m2 dzdw 



(27ri) 2 J To J Vl 2 e tz p m2 (z)w Xl+mi w(w - z) ' 



[z — l) m , m < rriQ 

(z-l) m °(z-2) m - m °, m>m . 

Moreover, Tq is a contour that encircles the pole w = but no other and 
is equipped with counterclockwise orientation. The contour ri 2 encircles 
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the poles z = 1,2 but no other and is also equipped with counterclockwise 
orientation. 



Remark 2.13. The expression (2.18) for the kernel can be obtained from 
the expression as given in [8] by the transform w i— > 1/w and zi-> 1/z. 

The fact that we deal with an determinantal process with an explicitly 
known kernel, allows us to analyze the system in detail. The main idea 
behind the proof of Theorem 2.10 is to write (h — E/i, <p) as a linear statistic, 



which allows us to rewrite the characteristic function at the left-hand side 



of (2.16) as a (Fredholm) determinant. The proof of Theorem 2.10 is then 
given in two steps: first we compute the limiting behavior of the variance 
of {h — E/i, 4>) and then we use the Fredholm determinant representation for 
the characteristic function to prove that the fluctuations are Gaussian. 
As we will show (see Proposition 4.3), the variance of (h —Mh, <p) can be 



expressed as a double sum involving the kernel K in (2.18). The asymptotic 
behavior for L — > oo can then be found by an asymptotic analysis for the 
double integral representation for K based on saddle point methods. How- 
ever, the disadvantage of this approach is that we also need to control the 
asymptotic behavior of the K near the boundary points, which then after- 
wards turn out to be negligible (note that the Gaussian free field is defined 
on a Sobolev space with Dirichlet boundary conditions). For this reason, we 
will follow a different approach that shows that the boundary is redundant 
in a more direct way. Instead of computing the asymptotics for K first, 
we provide an alternative expression for the double sum in the variance in 
terms of a quadruple integral and compute the (bulk) asymptotics after- 
wards. Moreover, based on the fact that we deal with Dirichlet boundary 
conditions, we give a probabilistic argument to show that the boundary has 
no effect on the fluctuations and that we can ignore the boundary in the 
Fredholm determinant identity for the characteristic function of (h — Kh, <f>). 
Apart from the conceptual benefit, this approach also reduces the amount 
of computational technicalities. Indeed, as the boundary has a complicated 
structure due to the different jump rates, giving case by case estimates on 
the kernel near the different parts of the boundary is cumbersome. 



2.6 Overview of the rest of the paper 

The rest of the paper is organized as follows. In Section [3] we prove Propo- 
sitions 2.4, 2.8 and 2.9 on the complex structure. In Section [4] we compute 
the limiting behavior of the variance of (h — E/i, <j>). In Section [5] we prove 
that the fluctuations are Gaussian and by combining this with the results 
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of Section [4] we obtain a proof for Theorem 2.10 Some of the statements 
in Sections [4] and [5] are based on steepest descent arguments that we post- 
pone to Sections [6] and [7j More precisely, Proposition 46 is proved in [7] and 
Lemma 5.3 is proved in Section [6} The proof of Theorem 2.7 is, although 
standard, given for completeness in Section [6j 



3 Proof of Propositions 2.4, 2.8 and 2.9 



In this section we prove Propositions 2.4, 2.8 and 2.9 



there exists a 



Proof of Propsition 2.4 ■ We prove that to every Q E 
unique (£, fi) E V such that Q = /j,) . Let us first consider (2.4) for 
fi < fio. By taking real and imaginary parts we can put F'(Q | £, fi) = in 
matrix form as 



(3.1) 



Re h Re 
Im n Im (n 




The determinant of this matrix can be easily computed to be 

Imft 



det 



Re n 



R e(n 



Im n Im n 




which is always strictly positive for E 
and solve (3.1 ) 

(3.2) 



1)| 2 ' 

Hence we can invert the matrix 



ni 



in - il 



in - il 



We recall that this formula is only valid for /j, < /xq. Concluding, for any 
n E HLl with |J7 — 1| 2 < fio/r there is a unique pair (£,/i) such that ft = 

n(6/x). 

By repeating the same procedure but now for \i > [j,q we obtain 



(3.3) 



2 



T + 



Mo + o I t 



MO 



Of course, this is only valid for fi > (jlq so that we need r- 



Mo 



|n-i|s 



> 0. Hence 



for every n E H_|_ with \Q — 1| 2 > fio/r there is a unique pair (£,//) such 
that n = n(£,//). 
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Concluding, we find that the map (£, fi) h-> fj,) is a one-to-one map 
from {(£,//) | /i < /io} onto 



|0 -IP < ^}, 

T 



and restricted to {(£, /x) E X> | fj, > Ho} it is one-to-one to 



l^-i| 2 >^}, 



Since these sets are disjoint and the union is H+, we see that the unrestricted 
map h-> fi"(£,/u) is indeed one-to-one from I? onto H + . 

The continuity is statement is immediate. □ 



Proof of Proposition 2.8 T he v ector {—dh/d^,—dh/d^i,\) is a normal to 
the graph of h. By Theorem 2.7 taking total derivatives and using F'(Q) = 
0, we obtain 



(3.4) 

and 

(3.5) 

where a 
have 



d - 1 d 1 

— h=-—lmF(n) = argO, 

— Im -F(O) = — arg(0 — a) arg f2, 

IT O/J, 7T 7T 



1 if /li < no and er = 2 if \i > /j,q. Hence by definition of 9j we 



(3.6) 



d 



-h = — 

<9£ 7T 

which proves the statement. 



and 



h n-o 2 



7T 



7T 



TT 



□ 



Proof of Proposition 2.9. By arguing as in the proof of Proposition 2.8 above, 
we obtain 



(3.7) 



d_ 



T 19 „, n . Im U 
h = -— ImF(fi) = . 

IT OT 7T 



sm &2 8111 (73 



1111 



(7 



By the law of sines we have 
(3.8) 

And hence we can write 

(3.9) ImO = |Q|sin0i =■ . 

o~ sm P3 a sm #3 

where a = 1 if /x < no and a = 2 if /x > fj,Q. The statement now follows by 
writing Im f2 and 9j in terms of the partial derivatives of h with respect to 
£, // and r using (13.41), (|3.5l) and (13.71). □ 



sin c/i sin (72 sin (71 srn( (/1 
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4 The variance of (h — E/i, 4>) 



The purpose of this section is to compute the limit of the variance of (h 
E/t, (A) as L — > oo. 



Proposition 4.1. Xe£ /i 6e as m (2.3), (j) aC 2 function with compact support 
in H + and (•, •) as in (2.14). 



Then 



(4.1) 



lim Var(/i 

L— >oo 



E/i, 



2 



where \\ ■ ||y is i/ie Sobolev norm ( |2.12 ). 



The proof will be given in Section 4.2. We will first provide a formula 
that expresses the variance of a general linear statistic in terms of a kernel R. 



We will show that (h — Kh, (j)) is in fact a linear statistic and Proposition 4.1 
then follows by the asymptotic behavior for R as L — > oo. The proof of the 
asymptotic behavior of R is based on a steepest descent argument and will 
be postponed to Section [7j 



4.1 The variance of a linear statistic 

Let / :ZxN— >Rbea function with finite support and define the random 
variable Xf as 



(4.2) 



Xf= /( x ' m )> 

(x,m)&C 



where C is a random configuration of points. 

Random variables of the type Xf are generally referred to as linear statis- 
tics and play an important role in the study of determinantal point processes. 
It is standard (and straightforward to derive from ( |2.17 )) that 



(4.3) 



and 



KXf = f(x,m)K(x,m,x,m), 

(x,m) 



(4.4) V&rXf = f(x,m) 2 K(x,m,x,m) 

(i,m)eZxN 

^ f(xi,m 1 )f(x2,rn2)K(x 1 ,m 1 ,X2,m 2 )K(x2,m2,x 1 ,m 1 ), 

(xi,mi)eZxN (z2,m.2)eZxN 



where K is the kernel in (2.18). We will rewrite (4.4) in terms of a kernel 



R. To this end, we need the following lemma. 
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m\ > 777-2 



mi < m2 



Figure 6: Deforming the contours such that the kernel K(xi, mi, x 2 , mi) 
is represented by one double contour only. The solid line represents the 
contour Tq. The dashed contour represents Ti^. The dots are the points 
0, 1 and 2 which are poles for the integrands. 



Lemma 4.2. With K as in ( |2. 18 ) we we have that 



(4.5) 



^ K(xi,m 1 ,X2,ni2)K(x2,m 2 ,xi,mi) 



K(xi,mi; x\, m\)8. 



mi,m2 



for (x\, mi) £ZxN and 777,2 G N. 



Proof. We start by rewriting the integral formulas for K in (2.18) as follows. 
If mi > ni2, we deform Tq such that it also goes around Ti^. Due to the 
term (w — z)" 1 we pick up a residue that results into single integral over 1^2 . 
However, by the assumption mi > m-2 the integrand of the single integral 
has no pole inside Tip- Hence it vanishes and we are left with the double 
integral only, where now Tq also goes around T^2- If on the other hand, 
mi < m 2 , then we deform the contour Ti 2 such that it also goes around Tq. 
The residue that we pick up in this way results into a single integral over Tq 



that exactly cancels the single integral that was already present in (2.18). 



Concluding we can rewrite (2.18) as 



(4.6) K(xi,mi,x 2 ,m 2 ) 



1 



(2^i) 2 



r Jt 



e tw p mi (w)z X2+m2 dzdw 
e tz p m2 (z)w xi+mi w(w — z) ' 



where Tq now also goes around I\2 if mi < m,2, and if on the other hand 
mi > m-2, we have that also goes around Tq. See also Figure [6j 

The next step in proving (4.5) is to write K(xi,mi,X2, m 2 ) and K(x2,ni2, x%, mi) 



both as one double integral. Note that we can always deform the contours 
such that all contours are circles around the origin and such that the radius 



of r' and Ti 7 2 are equal. After these preparations we substitute (2.18) into 
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the sum (4.5) and take the sum under the integral. The sum is then over 
terms -^Sr and we have 

W x 2 



(4.7) 



^ (2vri) 2 



e tw ' Pm2 (w') 



z 



X2+m,2 



dw'dz 



e tz Pm 2 (z) w 



X2+rri2 + l 



(w — z)(w' — z') 



2vri ./ r 



dz 



(w — z)(z — z') 



The integral at the right-hand side vanishes if w and z' are either both 
inside the contour Ti^ or both outside. These situations precisely occur 
when mi / m 2 - Hence we proved (4.5) for mi / m%. Finally, if mi = m.2, 
we have that z' is in the region enclosed by T12 and hence the left-hand side 
of (4.5), by using (4.7) and picking up the residue at z = z', equals 



(4.8) 



(2vri) 5 



p mi {w)z' Xl+mi dz'dw 



e* 2 p mi 1 



z )W 



w[w 



where To goes around Ti^'. Hence it equals K{x\,m\,x\,m\) by (4.6) and 
we also proved (4.5) for the remaining case mi = mi. □ 

Using the fact that the diagonal of K 2 and K agree, we can rewrite 
the variance of linear statistic in a useful different form as presented in the 
following proposition. 

Proposition 4.3. Define the difference operator D byDf(x,m) = f(x,m) — 
f{x - 1, m) for f : Z x N -»• R. Then 



(4.9) 



Df(xi,m 1 )Df(x 2 ,m2)R(x 1 ,mi,X2,m2), 



where 



(4.10) R(yi,mi,y2,m2) 

f Y,xi>m Yjx2<v2 K{x 1 ,m 1 ,X2,m2)K(x2,m 2 ,xi,m 1 ) yi>y 2 , 
\ T,xi<yiT l x 2 >y 2 K ( x ^^ m i^ x 2: m 2)K(x2,m 2 ,xi > rn 1 ) yi < y 2 . 

Proof. The right-hand side of ( |4.4[ ) consists of two terms. Let us first con- 
sider the first term. By writing f{xj,rrij) = Yl y < x Df and changing the 
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order of summation we get 



/l f(x 1 ,m 1 ) 2 K(x 1 ,m 1 ,xi,m 1 ) 



11,1711 



E ^2Df(y 1 ,m 1 )Df(y 2l m 1 ) ^ K(x 1 ,m 1 ,x 1 ,m 1 ) 

xi>max{yi,y2) 



2/1, mi t/2 



By separating the cases y\ < y 2 and y\ > y 2 and using (4.5) to introduce 
extra sums over x 2 and m 2 , we obtain 

(4.11) 

E f(xi,m 1 ) 2 K(x 1 ,m 1 ,x 1 ,m 1 ) = ^ ^ Df(y 1 ,m 1 )Df(y 2 ,m 2 ) 
xi,mi yi,miy2,m2 

Y.x 1 > yi Ex 2 K(x 1 ,m 1 ,x 2 ,m 2 )K(x 2 , m 2 , x 1 ,m 1 ), y x > y 2 
J2 X1 t2x 2 >y 2 K(x 1 ,m 1 ,x 2 ,m 2 )K(x 2 , m 2 , xi, mi), y 1 < y 2 . 

Note that after summing over x 2 , the summand vanishes if mi ^ m 2 . 
The second term in (4.4) is a double sum that we can rewrite as 

^ ^ f(x 1 ,m 1 )f(x 2 ,m 2 )K(xi,m 1 ,x 2 ,m 2 )K(x 2 ,m 2 ,x 1 ,m 1 ) 

= E E E D f(yi> m i) E D f(y^ m ^) 

xi,mi x 2 ,m 2 yi<xi y2<x 2 

x K(xi,mi,x 2 ,m 2 )K(x 2 ,m 2 ,xi,mi), 

which after changing the order of summation reduces to 

(4.12) ^ ^ f(xi,m 1 )f(x 2 ,m 2 )K(x 1 ,m 1 ,x 2 ,m 2 )K(x 2 ,m 2 ,x 1 ,m 1 ) 

E E D f(yii m i) D f(y2i m 2) 



xi,mi X2,m 2 



yi,mi V2,m2 

x ^ ^ K(xi,mi,x 2 ,m 2 )K(x 2 ,m 2 ,xi,m 1 ) 
xi>yi X2>y2 



Inserting (4.12) and (4.11) in (4.4) gives the statement. 



□ 



We end this paragraph with a symmetry property of R that is useful 
later on. 
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Lemma 4.4. If m\ ^ then 

(4.13) V] ^ K(xi,mi,x 2 ,m2)K(x2,m2,xi,mi) 

xi>yi X2<V2 

= ^2 ^2 K(xi,mi,X2,m2)K(x2,m2,xi,mi) 



xi<yi x 2 >V2 



Proof. The proof follows by applying (4.5) twice and changing the order of 
summation 

(4.14) ^2 E K(x 1 ,mi,x 2 ,m 2 }K(x 2 ,m 2 ,x 1 ,mi) 
xi>yi x 2 <y2 

= - E E K ( x i^ m ij x 2,m2)K(x2,m 2 ,xi,m 1 ) 

xi>yi X2>y2 

= ^2 E K ( x ij m ij x 2,m2)K(x2,m 2 ,x 1 ,m 1 ). 

Xl<yi X2>V2 

□ 



4.2 Proof of Proposition 4.1 



In this paragraph we prove Proposition 4.1 We need the following lemma 



that expresses (h, <f>) as a linear statistic. 

Lemma 4.5. Let h be as in (2.3), <f> a C 2 function with compact support in 



,, (h,(f>) as in (2.14) and f be defined by 

J2 Act>(n(y/L,m/L))J(y/L,m/L). 



(4.15) f(x,m) 



LJ 



y<x 

(y,m)eLT> 



Then (h, 0) = Xf. 

Proof. Let C be a random configuration of points. Then 



L 2 



I J2 h(y,m)Act>(n(y/L,m/L))J(y/L,m/L). 



(y,m)eLT> 



By inserting definition of h in (2.3) and changing the order of summation 
we get 



V- 



l E E &<t>(to(v/L,m/L))J(y/L,m/L), 



(x,m)GC y<x 

(y,m)£LV 



which is the statement. 



□ 
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Now that we have shown that (h, eft) is a linear statistic, we can express 



its variance in terms of R by Proposition 4.3 Hence it suffices to find the 
asymptotic behavior of R, that we present in the following proposition. 



Proposition 4.6. Fix 5 > and set 

(4.16) 



I Xj = [L£j] , 
\rrij = [Lfij], 

where [x] stands for the integer part of x. Then 



(4.17) 



lim R(xi,mi, X2,m,2) 

L— >oo 



2vr 2 



log 



0(£i,Mi) - 0(£ 2)/ u 2 ) 



o(£l,mi) - ^(6)^2) 



uniformly for fj,j) in compact subsets ofT> and Ml) - ($2> M2) || > L 2 
Moreover, we have 

(4.18) R(xi,mi,X2,m2) = 0(logL), as L — » 00, 

where the constant is uniform for (£j,Mj) in compact subsets ofD. 

The proof of Proposition |4.6| will be postponed to Section [7j Finally we 



prove Proposition 4.6 



Proof of Proposition 4-6 Write (h,4>) as Xf with 

f{x, m) = -^J2 A0(n(z/A m/L))J(x/L, m/L). 



y<x 



Then £>/ = -L~ 2 A^. Hence by Var(/i - E/i, 0) = Var(/i, 0) = Var X/ and 



by Proposition 4.3 we obtain 

Vax{h-Eh,<f>) = -j± Yl Yl &<P{tt(x 1 /L,m 1 /L))A(j){n(x2/L,m2/L)) 

(xi,mi) (x 2 ,m 2 ) 

x R(xi,mi,X2,m2)J(xi/L, mi/L) J(x 2 /L, m 2 /L). 



The right-hand side can be viewed as a Riemann sum and hence by the fact 
that 4> has compact support we have that 



lim Var(/i — E/i, 




X log 



A0(n(£i,Mi))A0(n(£2,02)) 

0(6, Ml) - 0(£ 2 ,^ 2 ) 



0(6) Mi) - 0(^2, M2) 



^"(6 , Mi ) J {£2 , M2 ) d£i d/ii d£ 2 d/U 2 . 
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Here we used (4.17) for points {£j,fJ>j) that are far apart and (4.18) to show 
that the contribution of the points that are close is negligible. By a change 
of variables we obtain 



lim Var(/i - Eh, > 

L— >oo 



A0(fia)A^(n 2 )— log 
Ztt 



fii - n 2 



fil - fi 2 



dm(J7i)dm(r22) ; 



where dm stands for the planar Lebesgue measure. The fact of the matter is 
that the logarithm in the integrand is the Green's function for the Laplace 
operator on the upper half plane with Dirichlet boundary conditions. There- 
fore 

<£(n)A0(n)dm(fi) = 



lim Var(/i — Eh, 

L— s>oo 



where the last step follows by integration by parts. 
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□ 



5 Gaussian fluctuations 



In this section we present a proof of Theorem 2.10. The starting point is 



again that we can write (h, 4>) as a linear statistic Xf with 

f(x, m) = -^J2 m/L))J(x/L, m/L), 



y<x 



as given in Lemma 4.5 Note that the support of / is unbounded (in con- 
trast to the compact support of <f>). Indeed, the function / is constant on 
horizontal rays that are at the right of the support of <ft (scaled with L). 

The unbounded support of / turns out to be inconvenient in the proof. 
Therefore, we split the function / = f\ + fa such that fa has bounded 
support containing the (scaled) support of Acft o £1. To this end, fix e > 
and define 



(5.1) 

Now split / as 
(5.2) 



D £ = {(^)GD:ImnK, /J )> £ }. 



/ = /l + /2 



fi = fxv e ■ 



where xv e stands for the characteristic function of the set T> £ . The point is 
to choose e small enough such that T> e contains the support of A^ o fi. See 
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Figure 7: The solid curve represents dT> and the dashed curve &D £ . The 
shaded region is the (compact) support of A0 o Q,. 



also Figure [7} In that case, the function fa contains all the information of 



In Section 5.1 we first prove that the contribution of Xf l to Xf is neg- 
ligible. More precisely, the variance tends to zero as e 4- after taking the 
limit L — > oo. Then in Section 5.2 we prove that the fluctuation of Xf 2 are 
Gaussian. Finally, based on these two results and a standard probability 



argument we prove Theorem 2.10 in Section 5.3 



5.1 The variance of X 



h 



For Xf x we prove the following lemma, for which the main ingredient is 



Proposition 4.6 



Lemma 5.1. For e > 0, let X^ and fa be as in (4.2) and (5.2). There 
exists a positive function g^, with lim £ |o 9<f>( e ) = 0, such that 



limsupVarXft < (g^(e)) . 

L— >OD 



Proof. By Proposition 4.3 we need to compute Dfa. We recall that <f> has 
compact support, so that we can choose e > small enough so that A0 o Q 
has a support that is entirely contained in V £ . In that case we have that 
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Dfi(x,m) 7^ only if (x,m) is in the set 



(5.3) B = {(x, m) £ LV £ \ (x - 1, m) G LB \ D £ } 

U{(x,m) E L(V\ V £ ) | (x - l,m) £ LV e }. 

Note that B is a discrete set of points close to LdT> £ . Moreover, 



(5.4) Dh{x,m) 



L 2 



A</>(Sl(y/L,m/L))J(x/L,m/L) 



y<x±l 



if (x,m) 6 B and Dfi(x,m) = otherwise. Hence 
(5.5) 

VarX fl = ^ ^ Df 1 (x 1 ,mi)Df 1 (x2,m 2 )R(xi,m 1 ,X2,m 2 ). 

(xi,mi)£B (x2,rri2)&B 

The right-hand side is a Riemann sum for a double integral over 9P £ x dT> £ . 

First note that the right-hand side of (4.17), for points (£j,/^j) € <92? e , 
can be written as 



(5.6) 



2vr 2 



log 



^(£1,^1) - ^(£2,^2) 



1 

1 



4vr 2 



log 
log 



1 + 



2ie 



1 + 



- ^(6,^2) 
4e 2 



|fi(6,/ii)-0(£ 2 ,M2) 



where in the last inequality we used that f2(£i,Mi) — ^(£21^2) S K for 

Mi) e 

In order to interpret (5.5) as a Riemann sum for an integral over &D £ x 
we use the following parametrization for T> £ that is based on the pro- 
jection onto the vertical coordinate. There exists a continuous function 



(e,t) H- {£e(t),(i e (t)) and points -00 = t < < t^' < t^' < i 4 
that depend continuously on e, such that fjf £ (t) is constant on (ty > ij-+x) and 
takes the values ±1, and 

2? e = {(&(*), M*)) |*GR}. 



By (5.4) we see that the limiting behavior of ±L_D/i([L£ e (i)], [L/i e (i)]) is 
given by 



(5.7) 



5(t) 



(*) 



A^(n(e,/i e (t))j(c,/i e (t))de. 
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By taking the limit L — > oo in (5.5) and using (5.6) and (5.7) we obtain 



lim sup VarX^ 

L— >oo 



< 



4tt 



mw&i log 



i+ 



4e^ 



- Q(s)| 2 



dsdt =: (g^e)) 2 - 



Note that since cj> has compact support, also 5 has compact support and 
hence the g^{s) < oo. By taking the limit, we have lim £ ± g^e) = and this 
proves the statement. □ 



5.2 Gaussian fluctuations for Xf 2 

The purpose of this paragraph is to prove the following proposition. 



Proposition 5.2. Let f 2 be as in (5.2). Then 

(5.8) hm (E (^A- H A)) - e -i* 2VarX /2 

uniformly for t in compact subsets of C 



The proof of Proposition 5.2 is based on some asymptotic results on 
the kernel K, that we will first present. We recall that for p > 1 the p-th 
Schatten norm of a matrix A is defined as 



(5.J 



= Tr \A*A\ p / 2 = J2 °j( A ) P , 



where Cj(A) are the singular values of A, counted according to multiplicity. 
Moreover, ||^4||oo = o~ max (A), i.e. the maximal singular value. 

Lemma 5.3. Let 5 > and D c be a compact subset of T> . There exists a 
function G : (Z x N) n LT> C — > (0, oo) such that the matrix 

(5.10) K G (xi,mi,X2,m 2 ) = —, — ' zK(xi, m 1 , x 2 , m 2 ) 

G(x 2 ,m 2 ) 

/or (xi, mi), (x2, m,2) G LP C satisfies 

1. \\K G \\ 2 = 0(L 1+S ) as L -> oo. 

g. ||if G || 4 = C(L)as L -> oo. 

5. H^gIIoo = O(L) as L -> oo. 
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Moreover, if we set £lj 
then 



Q(xj/L,rrij/L) and Fj(-) = F(- \ Xj/L,mj/L), 



4- for (xj,mj) £ LV C with ||(xi,mi) — (x2,m2 
(5.11) 
K(xi,mi,X2,m,2) - 



> L2 +<5 we have 



■(- 



e L(Fi(ni)-F 2 (n 2 )) 



2 7 iAn 1 (n 1 - fi 2 )L(-Ff(^ 1 )) 1 /2(^(^ 2 ))i/2 

e L(Fi(ni)-F 2 (n 2 )) 



+ 



^(flj - 2 )L(-F{'(fi 1 ))l/2(^(fi 2 ))l/2 

e L(Fi(fii)--F 2 (n 2 )) 
flip! - f] 2 )L(-F{'(n 1 ))V2(^(fi 2 ))i/2 



+ 



1 + o(l- 5 / 2 )), 



- 2 )£(-F{ / (Ol)) 1 /2(^(0 2 ))l/2, 

as L ^ 00, where the constant is uniform for (xj,mj) £ LT> C . The 
signs of the square roots (— ^"(fJi)) 1 / 2 and {F^f^l)) 1 / 2 are given in 
Lemma\6.1\, 



The proof of this lemma relies on a steepest descent analysis on the 
integral representation (2.18) for the kernel K. This analysis requires a 
significant amount of work and we therefore postpone it to Section [6} The 
choice in the sign of the square roots depends on a deformation of the contour 
in the steepest descent analysis. Since the the sign does not matter for the 



proofs in this section we omit it here and postpone it to Lemma 6.1 



In the proof of Proposition 5.2 we will also need the following. 
Lemma 5.4. Let f be as in Lemma (4.5) and f = /i + / 2 as in (5.2). Then 

1. II/2II00 = 0{L- 1 ) as L^ 00. 

2. [I/2H2 = 0(1) as L -> 00. 



Proof. This follows easily by the definitions of / 2 in (5.2) and / in (4.15), 
and the fact that II Aa 



< 00 by assumption. 



□ 



The last ingredient for the proof is the following result that is based on 
a lemma in 14 



Lemma 5.5. 

(5.12) 



We have that 



lim Tr(/ 2 iTr 

L— >oo 



> 3. 
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Proof. The proof of this statement follows by (5.11) and using the same 
arguments as in \14\ Sec. 7]. We will briefly highlight the main points here. 
To start with, expand the trace 

(5.13) TT(f 2 K) e = ••• f(xi,m 1 )---f(x e ,m e ) 

x K (xi, mi, x 2 , m 2 )K(x 2 , m 2 , x 3 , m 3 ) ■ ■ ■ K(xt, m e , xi, mi). 

Note that the support of f 2 is contained in LT> C for some compact subset of 
T>. Moreover, the contribution of points (xj, nij) that are close is negligible. 
To this end, we note that H/2II00 = 0(L _1 ) and that we can replace K with 
Kq and use (6.22), (6.25), (6.27) and ( 6.28| ) to show that the contribution 
coming from points for which ||( ) — (xj + i,rrij + i)\\ < L l / 2+& for some 
1 < j < ^, tends to zero as L — > 00. 

Hence in (5.13) we only need to consider points (xj / L,rrij / L) E D c for 



which the subsequent point are sufficiently far apart such that we can use 



(5.11 ). By expanding parenthesis we obtain 4 terms. Most terms are highly 



oscillating and therefore their contribution to the sum (5.13) is negligible. 



Hence there are only A£ surviving terms, that do not oscillate. Up to mul- 
tiplication with a symmetric function, each term can be written as 



(5.14) 



n 



1 



U)j+l 



where ojj = IX/ or Qj and w^+i = oj\. The point of the proof is that the sum 



(5.13) is invariant under permutation of variables. Hence we can replace Uj 



in (5.14) with (jJ c u\ for any permutation a. Moreover, we can replace it by 
a sum over any set of permutations. Lemma 7.3 in [TJ] tells us that the sum 
of (5.14) over ^-cycles is zero and hence we obtain the statement. □ 



We are now ready for the proof of Proposition 5.2 The proof relies on 



a Fredholm determinant identity for the characteristic function of a linear 
statistic. The statement then follows by estimates on the operator in the 
determinant. 



Proof of Proposition 5.2. First we note that since Xf 2 is a linear statistic, 
it is standard that we can write the characteristic function as a (Fredholm) 
determinant 



(5.15) 



E(expitX /2 ) = det 



(l + (e i</2 



l)K 
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Indeed, by writing the exponential of the sum as the product of the expo- 
nentials we obtain 



E(exp itX h 



E 



(x,m)eC 



where C at the right-hand side is a random configuration of points. By 
expanding the product we see that the right-hand side is a Fredholm deter- 
minant (see for example for more details). 

Note that fi has support in LT> C some compact D c CT>. Hence the same 
is true for e 1 '-^ 2 — 1. This implies that we can (and do) restrict the matrix 
K to the points on the grid Z x N that are inside LD C . In particular we can 



apply the results of Lemma 5.3 



Since the determinant is invariant under conjugation with any invertible 
matrix we have that 



(5.16) 



det (l + (e ii/2 - l)if) = det (1 + (e ii/2 - l)K G , 



where Kq is as in ( 5. 1Q[ ) . It is useful to replace Kq to K since we can control 
norms of Kq as L — > 00 as given in Lemma 5.3 In the remaining part of 
the proof, we analyze the asymptotic behavior of the determinant at the 



right-hand side of (5.16). 



Let us first make some remarks on the operator in the determinant. 



First, by Lemma 5.4 there exists a constant C\ such that 



(5.17) 



Jtf 2 



llloo < tC\L 



By combining this with the norms on Kq in Lemma 5.3 and using the fact 
that ||^4-B|| P < H^Hoo ||-B||p, we see that there exists constant C2 and C3 such 
that 



(5.18) ||(e it/2 -l)^ G ||oo<tC 2 and ||(e^ 2 



1)K G \U < tC 3 . 



Clearly, the constants Cj do not depend on L. 

Now we return to (5.16). The determinant can be rewritten as 



(5.19) 



det 1 + e 



exp 



E 

k£=1 



(-1) 



e+i 



Tr (e 



Jtfr 



1)K G 



This expansion is valid for t in a neighborhood of the origin in the complex 
plane. A priori, this neighborhood depends on L. Our first task is to show 
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that the neighborhood can be chosen independent of L. To this end we note 
that for £ > 4 we have 



Tr ((e 



Jtf 2 



l)K G 



< 



Jtf 2 



l)K G \\\ 



< 



Jtf 2 



m G \\t 



-4|| (e ith 



m G \\ 



where in the last inequality we iteratively used the inequality H^Hp+i < 



l^llpll^lloo valid for A in the p + 1 Schatten class. Now by (5.18) we then 



have that there exists constants C4, D > 0, independent of L, such that 



(5.20) 



Tr (e 



1)K G 



< fClD. 



This means that (5.19) is valid for \t\ < C 4 1 , for all L. 



The next step is to analyze the asymptotic behavior of the traces in the 



exponential at the right-hand side of (5.19). We claim that 
(5.21) lim Tr ((e itf2 - 1)K G Y = 0, I > 3. 

L— >oo V / 

To this end, we first observe that we have the following inequality 



(5.22) 



Tif^K G ---f^K G \ < \\f 2 \\^\\K G \\ 2 2 \\K G 



\l-2 
1 00 ' 



> 2. 



Here we used the well-known identities | Tr^4| < ||^4||i, ||^4i?||i < ||A||2||B||2 
and || Ai?||2 < ||j4||2||-B||oo where || • ||i stands for the trace norm. If s\ + S2 + 



• • • S£ > £ + 1 the trace converges to zero as L — > 00 by Lemma 5.3 l)+(iii) 



and Lemma 5.4 i). Hence, by expanding the exponential, we see that the 



only term in (5.21) that possibly remains in the limit, is therefore the term 
with Sj = 1. Hence we proved 



(5.23) 



lim ( Tr ( (e itf2 



1)K G ) -Tr(itf 2 K G )> 



0, 



£ > 3. 



Now (5.21) follows by (5.23) and (5.12). 



By (5.19) and (5.21) we have 



(5.24) lim det fl + (e^ 2 - 1)^ e --*(e«*-l)K G+ ±Tr((e«f2-l)K a ) = ^ 
L— >oo \ / 

uniformly for t in a neighborhood of the origin. Note that we also changed 



the order of the limits, which is allowed by (5.20). 
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We proceed by simplifying (5.24) a bit further. 

By the identities after (5.22) and Lemma 15.41 we have 



Trf i K G \<\\f\\ t - 1 U\\ 2 \\K G \\ 2 = 0(L 2 + s -% asL^oc, i > 3, 



and hence by expanding the exponential we obtain 



(5.25) 



5-U 



as L — > oo. Moreover, by (5.22) we have 



(5.26) Tr ( (e^ 2 - l)K G ) = -t 2 Ti{f 2 K G f + 0{L 



Inserting (5.25) and (5.26) in (5.24) leads to 



(5.27) 



lim det (l + (e^ 2 - l) K G ) ^ hK G+ \ Tr(flK G - h K Gh K G ) = L 



By (4.3) and (4.4) and the fact that we can replace Kq back to K in the 
trace we can write this as 

(5.28) lim det (l + (e itf2 - l) K G ) e ~ mx h+^ VarX ^ = 1, 

Since the variance of Xj 2 remains bounded as L — > oo (indeed, Xj 2 = 
Xf — Xj x and both VarXy and VarX^ are bounded), we also have 



(5.29) lim ( det ( 1 + ( e^ /2 -l)K G )e 



-HEX 



Si — e 2 



0, 



uniformly for t in a neighborhood of the origin. Combining (5.29) with 
(5.15) gives (5.8) in a neighborhood of the origin. 

To prove that (5.29) (and hence (5.8)) holds for all t G C we argue as 
follows. Since (5.25) and ( 5.26| ) hold for all t € C it is sufficient to prove 
that (5.24) holds for all i £ C. To this end, we note that 



(5.30) 

f L (t) := det (l + {e ith - l)K G ^j e-Ti(e ith - l)K G + - Tr ({e ith - l)K G 2 



det4 (l + (e «A _ 1) Kg } e kTr((e»h~l)K G ) 
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where det4 is the regularized determinant of order 4 (in fact, this is how the 
regularized determinant is defined [21] )■ For such determinants we have the 
inequality HQ Th. 9.2(b)] 

|det 4 (l + A)\ < expr 4 || J 4|| 4 , 

for some constant T 4 (independent of A). Applying this inequality with 
A = (e ith - l)K G in gjjg ) and using fl5T8] ) and we see that {f L } L is 

a sequence of entire functions that is locally bounded. By Montel's Theorem, 
it has a subsequence converging uniformly on compact subsets of C. By 



analyticity and (5.24), the limit of any convergent subsequence must be the 



constant function 1. Hence, the entire sequence converges to 1 and (5.29) 



holds uniformly for t in compact subsets of C. Combining this with (5.15) 
gives the statement. □ 



5.3 Proof of Theorem [2Tl0l 



Proof of Theorem 2.10, In view of Lemma 4.5 we can reformulate the state- 
ment as 



(5.31) 



lim Eexp(it(X/ - EXr)) = exp 

L— >oo 



V i 



t e 



with / as in Lemma 4.5 Split / = /1+/2 as in (5.2 ), define Yf j = Xf j —¥,Xf j 
and write 



(5.32) (Eexp^ty/) -exp(-^ 2 ||0|||7)| < |Eexp(itY» -Eexp(itY/ 2 
+ |E exp(ity /2 ) - exp(- \t 2 Var Y h ) \ 



+ |Eexp(-it 2 VarY /2 ) -exp(-ii 2 

We estimate the three terms at the right-had side separately. 

Since Yt and Yt 2 are real valued and ~KYr = EYr 2 = 0, we have 



01- 



(5.33) 



|Eexp(iiy / ) - Eexp(itY /2 )| < E |exp(itl7) - exp(itF /2 )| 
< \t\n\Y f - Y h \) < \t\Jv a r(Y f -Y f2 ) = |i|v/VarYft. 



Hence by Lemma 5.1 there exists a positive function gAe) with lim e ^o 9<t>{ e ) 
and 



(5.34) 



limsup |Eexp(itY/) — Eexp(itY) 2 )| < \t\g^,{e). 

L— >OD 



This estimates the first term at the right-hand side of (5.32) 
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Proposition 5.2 deals with the second term. 
As for the third term, note that 



VarF; = Var(y /l +Y h ) = Var Y /x + 2E(Y fl Y h ) + Var Y h . 



(5.35) 



Hence by applying Cauchy-Schwarz and Lemma |5.1| we see that there exists 
a constant A > such that 



(5.36) lim sup | Var Yf 2 — Var Yf \ 



L— »oo 



< limsup (2 v / Vary /l VarV/ + Var Y fl ) < Ag^e). 



Finally, inserting (5.34), (5.36), (4.1) and (5.8) in (5.32) leads to 



(5.37) limsup \Eexp(itY f ) - exp(-±ii<£||^)| < (A+ \t\) g< p(e), 

L— >oo 

Since this holds for arbitrary e > we can take the limit e 4 and we arrive 
at the statement. □ 



6 Asymptotic analysis of K 



In this section we prove Lemma |5.3| and Theorem 2.7 The proofs are based 
on classical steepest descent arguments on the integral representation for 



the kernel (2.18), which by (2.4) we can write as 



K(xi,mi,X2,m2) 



(6.1) 



1 
2vri 



^L(F 1 (w)-F 2 (w)) 



dw 



(27Ti) 



MF 1 (w)-F 2 (z)) dzdw 
w(w — z) 



Here we used lower index Fj to denote the function F in (2.4) with parame 



ters (fij/fJ-j) = (xj/L, rrij/L). We will also write Qj to denote £l{xj/L } rrij/L). 
We will only need the asymptotics for (x, m) G LT> C for a compact subset 



T> c C T> in the proof of Lemma 5.3 For more background on steepest descent 



method applied to double integrals, we refer to [19 



6.1 Saddle points of F 

We start by investigating the saddle points of F, as well as the paths of 
steepest descent and ascent for Re F leaving from these points. 
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Figure 8: Contours of steepest descent/ascent for HeF leaving from the 
saddle points for the situation fi > //q- The black squares represent the 
points w = 0, w = 1 and w = 2. The black dots are the saddle points. 

If fj, < fj,Q then F has either two real saddle points or two complex 
conjugates and no more. If fi > fio there is an additional real saddle point 
in the open interval (1,2). We recall that T> is defined as the set of all pairs 
(£, fj,) such that we have complex conjugate solutions. Here we are interested 
in asymptotics for K for points (xj,m,j) £ LD C and hence we restrict our 
attention to the situation where there are two complex conjugate saddle 
points (and one real saddle point in the interval (1,2) in case fi > 

Since F is analytic in the upper half plane, the paths of steepest descent 
for Re.F are (parts of) level lines for ImF. The non-real saddle points 
are simple, resulting in a quadratic behavior for F near these points. This 
implies that we have two paths of steepest descent leaving from each non- 
real saddle point. One of the paths ends up in either w = 1 or w = 2 and 
the other ends in oo and is asymptotically parallel with the negative part 
of the real axis. There are also two paths of steepest ascent leaving from 
the saddle points. One of them is ending in w = and the other in oo 
asymptotically parallel to the positive part of the real axis. 

Finally, if /i > /iq there is a third saddle point in the interval (1,2). The 
paths of steepest descent are just the parts of the interval (1, 2) to the left 
and right of that point. The paths of steepest ascent end either in w = 
or oo. If the paths of steepest descent of the non-real saddle points end in 
w = 1 then the path of steepest ascent leaving from the saddle point in (1,2) 
ends in oo. Otherwise it ends in w = 0. See also Figure [8j 

6.2 Deforming the contours 

The next step in the steepest descent analysis is to deform Tq and to be 
paths of steepest descent /ascent for Rei 7 leaving from saddle points Q\ and 
r?2- However, by deforming the contours in this way we have to take the 
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term (w — z)~ l into account. Indeed, the deformed contours Tq and Ti^ will 
intersect at a point, say r/, in the upper half plane (and hence by symmetry 
also at fj in the lower half plane) . See also Figure [9j This implies that after 
deforming the contours we pick up a residue and we obtain 

K(x 1 ,m 1 ,x 2 ,m 2 ) = — [ e ^( w )-F 2 M)dw 

27T1 J(j W 

I eM*i(*)--F!.M). dzdw 



(2vri) 2 J ro J Tl2 w(w - z)' 

Here is a path connecting fj with r/ that crosses the real axis on the 
positive side if m\ > m 2 and on the negative side if m\ < m 2 . In both 
cases, the orientation of the path is from fj to rj. 
Let us introduce the notations 

(6.3) h = —( e L ^~ F ^ — , 

2vri J Cri w 

and 

(6.4) J 2 = J_/ / e L(F l{w) -F 2 ( Z )) ^dw 



(2vri) J ro J ri 2 w(w - z) 

and deal with the terms I\ and I 2 separately. 

6.3 Analysis of the double integral I 2 

First we consider the integral I 2 . In the first lemma that we now present we 
assume that f^i and Q 2 are not close to each other. 



Lemma 6.1. Let 5 > and T> c C T> compact. If (xj,rrij) € LT> C and 

, then 

I , e L(^i(fii)-i?2(n 2 )) 



-+s 

(xi,mi) — {x 2 ,m 2 )\\ > L2 , then 



(6.5) I 2 (x 1 ,m 1 ,x 2 ,m 2 ) 



+ 



2tt Vn^fi! - fi 2 )L(-F{'(0 1 ))V2(^(^ 2 ))V2 



ni(fii - 2 )L(-if (fii))i/2(f*'(n 2 ))i/2 

+ 



e L(F 1 (n 1 )-.F 2 (Q 2 )) 



fiiCHi - n 2 )L(-i^'(S 1 ))i/2(i^'(n 2 ))i/2 

+ - Ul + 0(L- S/2 )) 

^n 1 (U 1 -n 2 )L(-Fi'(n l ))^(F^n 2 ))^J { ^ { ,h 
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GO - <Xj + ' 

J-0 1 1,2 ^ 

Figure 9: When deforming the contours To and Tip in the double integral 



in (6.1) we pick up a residue that leads to an additional single integral. 



as L — )• oo. The constant in O is uniform for G T> c and does not 

depend on 5. The square root (—F^Qi)) 1 ^ 2 is chosen such that the orienta- 
tion of the line w = f^i + -, — „,j! u1 / 2 x tangent to Tq at Q.2, when traversed 
from —oo to oo, coincides with the orientation of Tq at f2i. Similarly for 

The square root (.Fg (f^)) 1 ' 2 is chosen such that the orientation of the 
line w = Q.i + (jr"^)) 1 / 2 ) ^ an 9 en ^ t° ^1,2, when traversed from —oo to oo, 
coincides with the orientation of Tip at Q.2- Similarly for {Fl^ (f^)) 1 ^ 2 - 

Proof. The point of the steepest descent method is that the main contri- 
bution of the integral comes from small neighborhoods around the saddle 
points. The parts of the contours To and Tip that are sufficiently far from 
the saddle points only give a contribution that is exponentially small. To 
obtain the leading order terms, we introduce local variables around the sad- 
dle points. Since we have a double integral and each integral gives rise two 
(conjugate) saddle points, there are four of such leading terms. 

Let us consider one of these four terms, namely the parts of the integrals 
around the saddle points fii and fl2- In small neighborhoods around these 
saddle points we slightly deform the contours To and Tip to their tangent 
lines and introduce the local variables 

(6.6) w = Qi -\ — -j= — and z = Q.i + 



£(-if(J2i))i/2 Vl(F^ 2 ))V2' 

for t G [-L 6 / 2 , L 5 / 2 ] and s G [-L 5/2 , L 5 / 2 ] where 5 > is the fixed constant 
in the Lemma. Note that with these new variables we have 



(6.7) 



lim L(Fi(w{t)) - Fi(fii)) = -\ and lim (F 2 (*(*)) - F 2 (Q 2 )) 

L— >oo L— >oo 



uniformly for t G [— L s / 2 , L s / 2 ] and by continuity also uniformly for fij) G 
V c . Note that at the endpoints t = ±L S we have that e L(Fl(, " W) " Fl(ni)) is 
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exponentially small as L — > oo. Moreover, from these points the contours 
are continued to paths of steep descent /ascent. This implies that the parts 
of the double integral that are not in the neighborhood of tt\ and il 2 (and 
their conjugates) is exponentially small as L — > 00. It therefore remains to 
compute the asymptotic behavior of the integral 



(6.8) 



(2m)*L(-F[' (fi 1 ))V2(^'(0 2 ))i/2 J_ LS J_ L s 



MFi(™(t))-F2(z(s))) _ 



dsdt 



w(t)(w(t) - z(s)) 



The main issue that remains is to handle the term l/(w(t) — z(s)). We claim 
that since Q± and Q 2 are sufficiently far apart by assumption, we have 



(6.9) 



1 



w(t) — z(s) fll — VL 2 



(1 + 0(L- S / 2 )), L^oo. 



in the new variables. Indeed, from (3.2) and (3.3) it is not hard to check 
that the inverse map 1— > (£, fi) is Lipschitz and hence there is a constant 
A such that 



(6.10) 



Combining this with the fact that F"(Q) is bounded from below for (£, fi) E 
T> c we have 

- ^ 2 ||F" 2 (fii i2 )| 1/2 \/Z > AL 5 



for some constant A. Hence 



\n l -n 2 \\F{' 2 {n l}2 )\ l / 2 VL 



< A-l L S/2 



for t G [— L 5 / 2 , L* 5 / 2 ]. This proves (6.9). Note that by continuity the constant 
A^ 1 can be chosen independent of (£j,[J<j) G F) c and hence depends only on 
V r but not on 5 > 0. 



Now substituting (6.9) into (6.8) leads to the following 



1 



(2vri) 2 L(-^'(0 1 ))i/2(^(r 22 ))i/2 J_ L& J_ L& 

e i(Fi(n x )-F 2 (n 2 )) 



,L{F 1 (w(t))-F 2 (z{s))) _ 



dsdt 



w(t)(w(t) - z(s)) 



(27ri) 2 L(-F^(n 1 )) 1 / 2 {F^(n 2 )) 1 / 2 {n 1 - n 2 )nx 



-LS J-L* 



MF 1 (w(t))-F 1 (n 1 ))-(F 2 (z(s))-F 2 (n 2 ))) dsdt ^ 1 + Q( L S/2s\ 
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By taking L large and computing the Gaussian integrals gives 

1 [ L [ L e L(F 1 ( W (t))-F 2 (z(s))) ^dt 

(27ri)2L(-^(I2 1 ))V2(^(n 2 ))V2 J_ L s J_ L S w(t)(w(t) - z(s)) 

e L(Fi(ni)-F 2 (n 2 )) 

-{l + 0{L- & / 2 )). 



2ttL(- j p{'(Oi)) 1 /2(^(o 2 ))i/2(o 1 - n 2 )ni 

The three other parts coming from the conjugate points Q\ and f2 2 can 
be dealt with in the same way. This four terms together give the leading 
asymptotic term as given in the statement. □ 

If Q.\ and VL 2 are sufficiently close, we need the inequality that is formu- 
lated in the following lemma. 

Lemma 6.2. Let 6 > and V c C V compact. Then there exists a constant 
such that 

(6.11) \h(xi,mi,x 2 ,m 2 )\ < • 

1,0 

for (xj,mj) 6 LV C with ||(xi,mi) — (x 2 ,m 2 )|| < L2 



Proof. The proof follows by the same approach as in Lemma 6.1 However, 



in this case (6.9) is no longer valid and we need to estimate this term in a 
different way. Note that since are close, we also have that F['(Qi) = 

F%(n 2 )(l + 0(£~ 1/2+<5 )), where the constant is independent of (Cj^j) G £>c 
as long as they are close. But then (— iq(Sli)) 1 / 2 = ±i(i ? 2 (f2 2 )) 1 / 2 where the 
±-sign depends on the choice of the square roots. Hence we have 

ii i 

(6.12) 



w(t) - z(s) v / L(-F 1 (0 1 )) 1 /2 n 1 -Q 2 + t±is + 0(L-V2+<5) • 

This means the double integral has a possible singularity (note if (£i,yUi) = 
(£2,^2) it certainly does), however this singularity is integrable. By pro- 
ceeding as in Lemma |6,1[ computing the Gaussian integrals and using the 
fact that F"(Q(£,{i)) is bounded from below for (£,//) E V c , we arrive at 
the statement. □ 

6.4 Analysis of the single integral I\ 

Next we give a bound for the single integral. 
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Lemma 6.3. Let I\ as in (6.3) and T> c <ZT> compact. Then there exists a 
constant C > such that 



(6.13) 



\Ii(xi,m 1 ,x 2 ,m 2 )\ < 



1 + R 



for (xj,mj) 6 LV C , where R = (xi — X2) 2 + (mi — m2) 2 ■ 

Proof. Let us first consider the case mi > m 2 . In that case, deform the 
contour to an arc from 77 and 77 that is a part of the circle centered around 
the origin with radius \r\\. A short calculation using (2.4) and the fact that 
mi > m,2 shows that 



(6.14) 



d. 
dt 



Re (Fi(r/e _i *) - F 2 (r]e~ [t )) < 0, < t < argr/. 



It also follows by the geometric principle that the distance of the point r/e -1 * 
to w = 1 and w = 2 decreases when t increases from to aigrj and attains 
its minimum at t = argry. It implies that 



(6.15) 



1 

2tt 



MFi(w)-F 2 (w)) dw 



W 



< LRe(Fi(r,)-F2(r,))_ 



We will refine the inequality for points (xj,nij) that are far apart using 
standard steepest descent arguments. Set 

(xi — X2,mi — 777.2) = -R(cos </>, sin 4>) 

and assume R > 0. Then write 

(6.16) L(F x {w) - F 2 {w)) = RGi, 2 (w, <f>). 

so that 



(6.17) 



1 

2^ 



C n 



dw 



w 



1 

2^ 



3 .RGi, 2 (uv 



I G?U7 



Note that 77 depends continuously on (^j,fJ>j) € X> c - By compactness of 
D c , there exists an r > be sufficiently small so that Im(7ye _ir ) > for all 
(xj, rrij) G LT> C . By ( 6.14[ ) there is a constant A > such that 



(6.18) 



Re (G^afae-* <f>) - £1,2(77)) < -At, < t < r. 



Again by continuity and compactness, the constant A can be chosen inde- 
pendent of <j> and 77, depending on P c only. 
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Finally, we split the arc C„ in three pieces, two small parts near rj and 



rj and a third remaining part. A simple estimate on (6.17) gives 

, dw 



(6.19) 



1 




2^ 




1 


,. r 


< — 




~ 2?r 


Jo 



3 HReGi, 2 (wW 



IV 



gRReGi^e-"^) + e RG ly2 (ve it ,<f>)\ df + ^-ArR^LKeG^ri 



For the second term in the second term at the righthand side we used (6.14), 



(6.16) and (6.18). To analyze the two parts near rj and 77 we use (6.18) and 



obtain 
(6.20) 



e flReG(r,e-«) dt < e RReG(rj) ( , 

Jo 



RAt 



dt < 



3 BReGi, 2 (rv 

AR 



Substituting (6.20) in the right-hand side of (6.19), inserting (6.16) and 



combing the result with (6.15) gives that statement for the case mi > rn-2. 

If mi < m2, then the statement follows by a similar argument. The 
main difference is that we deform C v to the arc of the circle with radius \r]\ 
that lies at the left of the origin. □ 



6.5 Proof of Lemma 15.31 

Proof. The function G that we use to define Kq as in (5.10) is given by 

3 LReF(f!) 



(6.21) 



G(x, m) 



Let us first give some estimates on Kq that can be deduced from the esti- 
mates on K that we have obtained so far. To this end we split 



(6.22) 
where 

(6.23) 



K G 



h,G + h,G 



T i \ G(x 2 ,m 2 ) , 
I j: G{x 1 ,m 1 ,x 2 ,m2) = mi, x 2 , m 2 ), 

Lr[Xl,mi) 



where Ij are as defined in (6.3) and (6.4) 
o 

(6.24) |/ 1>G | < 



From Lemma 16.31 we deduce 
A 



LRe(F 1 (t 1 )-F 1 {ni))-LRe(F2(v)-F2(n 2 )) 



1 + R 
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where R = ||(xi,mi) — (a?2j^2)||- Note that since r\ is an intersection point 
of Tq and ri j2 which are paths of steepest descent/ascent we have that 

Re (Fi(r7) - F 1 (Qi)) - Re (F 2 (ry) - F 2 (Q 2 )) < 0. 



Hence, from (6.24) we in have in particular 

A 



(6.25) 



\h,a\ < 



1 + R 



Suppose that R > L x / 2+5 . Then also and fl 2 are far apart and hence rj 
is far from either f^i or Q 2 . More precisely, there exists a constant B such 
that I r? - Oj| > BL- 1 / 2 ^ for j = 1 or j = 2 (see also (|6.10[ )). From ( |6.24[ ) 
and the quadratic behavior of -F near the saddle point, we obtain that there 
exists a constant B > such that 



.4 



1 + R 



(6.26) |/ liG | < 

for (xj,nij) G LP C . 

As for /2,G we deduce from Lemma's 6.1| and |6.2| that 

( 6 - 27 ) \l2,G(x 1 ,m 1 ,x 2 ,m 2 )\ < r u ^ . . 

— 12 2 | 

for (xj,rrij) G ££> c with ||(xi,mi) — (^2,^2)11 > ^ 1 ^ 2+<5 - Moreover, 



-4-1 



(6.28) 



\l2,G(xi,mi,x 2 ,m 2 )\ < 



A 2 



'L 

for (xj,rrij) G LP C with ||(xi,mi) - (x2,m 2 )|| < L 1/2+5 . 
Now we come to the proofs of (i), (ii), (hi) and (iv). 

(i). We split the Hilbert-Schmidt norm into two parts using the trian- 
gular inequality 1 1 1 1 2 < ll^i.db + 11-^2, ch- 

Let us start with the term g|| 2 - Then by (6.25 ) and a short calculation 
we obtain 



(6.29) ^2 \h,G(xi,m 1 ,x 2 ,m 2 ) 

(xi,mi)&LT> c {x2,m2)(zL'D c 



0(L 2 logL), 



as L — > 00. 
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Next we analyze ||/2,g||2- If (xi,mi) and (x 2 , rn 2 ) are close, then we have 

\I 2>G (x l ,m 1 ,x 2 )\ 2 = 0(L 2+s ), 



by (p8) 



(6.30) Yl 

(xi,mi)£LV. 



E 



(x 2 ,m 2 )£LT> c 

(a;i,m 1 )-(x 2 ,m2)||<L 1 /2+<S 



as L — > oo. 

Finally, we consider points (xi,mi) and (x 2 ,m 2 ) that are sufficiently far 
from each other. Then we use (6.27) to obtain 

(6.31) 



E 



E 

(xi,mi)eLT> c (x 2 ,m. 2 )eLT> c 

\\{x 1 ,m 1 )-{x 2) m 2 )\\>L 1 / 2 + s 



* E 

{x\,m{)(zL'D c 



h,G{xi,rni,X2,m 2 )[ 



E 



(i 2 ,m 2 )E iI 'c 

(x 1 ,m 1 )-(x 2 ,m 2 )\\>L 1 / 2 + s 



L 2 |^i - fi 2 | 



We interpret the right-hand side as a Riemann sum and estimate it by the 
integral 



(6.32) Yl 

(xi,mi)eZ/D c 



E 



(x 2 ,m 2 )£LT> c 

\(x 1 ,m 1 )-(x 2 ,m 2 )\\>L 1 / 2 + s 



< LA 




n(v c )x(n(v c )\B ni ) 



\h,G(xi,m 1 ,X2,m 2 )\ 
1 

-dm(Qi)dm(Q 2 ), 



|Oi-fi 2 | 2 



where Bq 1 is a ball around with a radius that is of order L 1 / 2+<5 , A > 
is some constant and dm stands for the two dimensional Lebesgue measure. 
Note that we also used that, by continuity and the fact that the Jacobian 
does not vanish, the Jacobian in (2.15) is bounded from below and above 
on compact subsets. As L — > oo the integral at the right-hand side grows 
logarithmically, so in particular we have 

(6.33) Y \h,G(xi,m 1 ,x 2 ,m 2 )\ 2 = 0(L 2 logL), 

(rri,mi)eLZ5 c (x 2 ,m 2 )&L(T> c \Bci 1 ) 



as L — > oo. The statement now follows by combining (6.29), (6.30) and 
(16331). 
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(ii) The proof is similar to the proof of (i). Write 



(6.34) 



\K g \\\ = TtK g K g K g K g 

, m2,xi,mi)K G (x2,m2, x 3 , m 3 ) 
x K G (x 4 , m 4 , x 3 ,m 3 )K G (x4, m 4 , x 1 ,m 1 )\, 



and insert ( 6.25 )-( 6.28). A short computation shows that the contribution 



of points that are close are negligible (in contrast to the double sum in the 
Hilbert-Schmidt norm) and the leading term comes from points (x 



j- 



m ; 



that are far apart. For such points we use (6.26) to deduce that the term 



Ii g only leads to an exponentially small contribution. Hence by (6.27) and 
arguing as in (i) we see that there exists a constant A > such that 



\K G \\i < AL 4 



dm(Qi)dm(fl2)<im(Q 3 )dm(0,4 

^1 — ^2||^2 — ^3||^3 — 1 1 _ 



Despite the singularities in the integrand, the latter integral is bounded and 
we arrive at the statement. Note the difference with (i) where the integral, 
due to the singularity, grows logarithmically with L. 



(iii) This follows from (ii) and the fact that 



> 



for p > 1. 



iv) In Lemma 6.1 we proved the statement for I2 instead of K. However, 



by adding the term I\ into the right-hand side of (5.11) and observing that 
I\ G is exponentially small for points far apart, we see that the statement is 
also true for K. □ 



6.6 Proof of Theorem I27T1 



Proof of Theorem 2.7. The proof is fairly standard for determinantal pro- 
cess with a kernel that is represented by a double contour integral. We will 
therefore allow ourselves to be brief. 

Write K = I\ + I2 with L as in (6.3) and (6.4). On the diagonal the 



main contribution comes from 1% instead of fy- Indeed, if (£, //) G T>, then by 
(6.28) and the fact the I^^g = h on the diagonal, we see that /2(C) M> £> A 4 ) 

On the other hand, from (6.3) and the fact that 



tends to as L — > 00. 
i] = fj) we have 



(6.35) 



2vri 



n dw 
n w 



argft(f,/i) 



7T 
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Figure 10: Configuration of contours in the quadruple integral representa- 
tion (7.1) for R in case m\ < mi- 



Hence we proved, with (x,m) as in (2.7), that 



(6.36) 



lim K(x,m,x,m) 

L— >oo 



7T 



for (x, m) such that (£,//) G P. 

From the proofs as presented in this paper, we only have (6.35 ) uniformly 
for (£, n) in compact subsets of V. However, by a steepest descent analysis 
for the points that are near the boundary dT> or outside T> it can be shown 
that (6.35) holds on T> and, if (£,//) is point in the upper half plane to the 
right of T>, then K(x, m, x, m) — > uniformly as L — > oo. 

Since 

i 

—Kh(y, m) = — \^ K(x, m, x, m), 

/ j ' > 



x>y 



we obtain from (6.36) 



lim —Mh(y,m) 



arg(0(e / ,/i))de / 



By using (3.4) this integral is easily computed and we obtain the statement. 

□ 



7 Asymptotic analysis of R 

In this section we prove Proposition [4T6J First we show that the kernel R can 
be expressed as a quadruple integral. The asymptotics for R then follows 
from steepest descent arguments on this quadruple integral representation. 
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Proposition 7.1. With R as defined in (4.10) and p m as defined in (2.19), 
we have that 

I f f f f e t ( w + w ') 

(7.1) R(yi,mi,y 2 ,m 2 ) = 



Pm»Pm 2 (w') z S/2+m2^i+mi dz'dw'dzdw 
Pm 2 (z)Pm 1 ( z ') W'y2 +rn 2wyi +m i (W - z)(w - z')(w' - z)(w' - z')' 

If mi < mz, then the contour Tq goes around w = 0, the contours T12 
and T' 12 9° around z,z' = 1,2 and Tq, and finally, the contour T' Q goes 
around w' = and T[ 2 and Ti )2 . 

If mi > 771,2, then the contour T' Q goes around w' = 0, the contours T^ 
and r' 12 go around z,z' = 1,2 and T' , and finally, the contour Tq goes 
around w = and T'± 2 and ^12- 

If mi = 777,2, then the contours T12 and T' 12 9° around z,z' = 1,2, the 
contours Tq and T' Q go around the origin and Fx 2 and T' t 2 respectively. In 
addition, if yi < y 2 , then T[ 2 goes around Tq and if, on the other hand, 
yi > 2/2; then Ti t 2 goes around T' . 

All contours have counterclockwise orientation. 

Proof. First we recall that, by deforming the contours To or r'12, we can 
write the kernel K as one double integral (and no single integral). Indeed, 



by arguing as in Lemma 4.2, we deform I 1 ^ such that it also encircles To 
in case mi < 7712 and, in case mi > 7712, we deform To such that it encircles 
1^2- See also Figure [6j This means that we can write the product of the 
kernels clS EL quadruple integral 

(7.2) K(xi, mi, X2, 7772 )K{X2, m2, Xl,X2] 



(2vri) 4 Jp J Vl 2 J r > o / r / , 

S W+W '^ p mi (w)p m2 (w') z *2+m2 z >x 1 +m 1 dz'dw'dzdw 

e t(z+z>) Pm2 (z)p mi (z') ro tt2+m2+i ro xi+mi+i ^ _ ^J7y _ z /y 

The location of the contours depends on wether mi < 7772, 777,1 = m2 or 
mi > m 2- Let us first assume that mi < 7772. Then the contour To goes 
around the pole w = 0, the contour T\ t 2 goes around z = 1 and z = 2 and 
Tq, the contour T[ 2 goes around z 1 = 1 and z' = 2 and r' goes around 
w' = and T\ 2 . 



There are two ways to compute R. By (4.10), (4.13) and the fact that 



mi / m 2 , we can either compute E Xl > yi E X2<y2 or £ xi<lfl T, X2 > y2 over 



the terms in (7.2). We choose to use the second way of summation. In that 
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Figure 1 1 : Possible configuration of steepest descent / ascent paths for Re Fj 
leaving from the saddle points Qj (dots). 



case, we make sure that \w'\ > \z\ for w S r' and z E and \w\ < 



for w G To and z' £ r' 12 . Inserting (7.2) into (4.10), changing the order of 
summation and integration, and using 



(7.3) 



1 x z' Xl 1 , 1 

— > = — , and — > 

w w xi z' — w w 



y X 2 



1 



xi <yi 



W * — ' W 

X2>V2 



1X2 



ur 



gives the statement in case m\ < m2, with the location of the contours as 



in Figure 10 



The situation m\ > 771,2 follows by the same arguments or by the sym- 



metry in (4.13). In case m\ = 1712, we do not have (4.13). As a result we can 



not choose between the two different ways of summing (7.2). The precise 



way of summing is given (4.10) and depends on wether y\ < 1/2 or y\ > yi- 



The statement then follows by the same reasoning as above for the case 

777i < m 2- D 

Remark 7.2. The locations of the contours in Proposition are chosen for 



convenience in the proof of Proposition 4.6 There are other deformations 



of the contours possible such that (7.1) holds. 



We are now ready to prove Proposition 4.6 



Proof of Proposition 4-6 The proof of (4.17) is an elaborate exercise in 



steepest descent arguments on the quadruple integral representation in (7.1), 
using the same type of arguments as in Section |6j 

Let us first present the main strategy. As in Section [6j we deform the 
contours of integration to paths of steep(est) descent/ascent for Rei^- leav- 
ing from VLj and Vtj (for a discussion on these paths see Section [6jl). Each 
time we deform one of the contours, we possibly pick up a residue due to 
the intersection with another contour. By repeating this a number of times 
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Figure 12: The paths of integration for the two double integrals after de- 
forming the contours. 



we are left with several integrals, each of which is an integral over paths of 
steep(est) descent /ascent of the integrands. Whenever the integrand con- 
tains an exponential, the integral tends to as L — > oo. There is only 
one integral that does not contain an exponential and this gives the main 



contribution and the right hand side of (4.17). 

Now let us discuss this procedure in more detail for the case mi < m-i. 
The case m\ > follows by precisely the same arguments but with 
replaced by T'^ and vice versa. In case mi = ni2, the initial contours for R 
are slightly different, but this does not lead to essentially new complications. 



We set mi < m-2 and start with the contours as in Figure 10 Let To 
first be a very small contour around the origin and T' Q a very large contour 
close to infinity far away from the other contours. Deform the contours T^ 2 
and Ti 2 to paths of steep (not necessarily steepest) ascent leaving from Q± 
and O2 and their conjugates, such that they do not intersect with Tq and 
r' (yet). Then deform the contours To and T' to paths of steep descent 
leaving from Cli and ^2 and their conjugates. By doing so, we pick up 
several residues, and are left with quadruple, triple and double integrals. In 
some of these integrals, we need additional deformations, but first we will 
collect the various terms that we obtained. 

For illustration purposes, let us also assume that for the deformed con- 
tours we still have that T' goes around Tq and that Tip goes around T[ 2 
(the other situations follow by similar arguments). See also the left picture 



in Figure 11, where Q2 is the top saddle point. In this case we are left with 



seven integrals: one quadruple integral (the same as in (7.1) but now over 



paths of steep ascent /descent), four triple integrals over contours that are 



indicated in Figure [13^ and two double integrals given in Figure 12 



We discuss the seven integrals that are obtained in this way, starting 



with the double integrals in Figure 12 The right picture represents the 
integral 

-1 pQi pQ 2 1 

(7.4) ^ / / j^dzdz'. 

V ' (27ri) 2 7 ni k (z-z' f 



45 




Figure 13: The contours for the four triple integrals after deforming the 
contours. 



The integral is easily computed and gives the Green's function at the right- 



hand side of (4.17). The double integral over the contours in left picture of 



Figure 12 is given by 



(7.5) 



(2vri) 



\\2 



m rm p l(f 1 (z)+f 2 (z')) 



Vl J V 2 



_^ e L(F 2 (z)+F 1 (z')) ( z _ z iy 



dz'dz, 



where rji and r/ 2 are the point of intersection in the upper half plane of 
ro and Fi t 2, and the contours T' and T' 12 respectively. Deform the paths 
of integration to arcs of circles centered at the origin with radius \rjj\. As 



shown in the proof of Lemma 6.3, these arcs are paths of steep descent and 
ascent for Re(i ? i — F2) leaving from 771 and r\2 (here we use the fact that 
mi < m 2 ). Moreover, since the m are intersection points of contours of 
steep descent / ascent for Re Fj , we also have 



(7.6) 

Re (Fifai) 



F 2 (r 1l )) < Re (Fi(fii) - F 2 (fi 2 )) < Re (Fi(?? 2 ) - F 2 (r ] 2)) . 



As a result, we see that the integral in ( 7.5 ) is exponentially small as L — > 00, 
when CI 1 and f2 2 are sufficiently far. 

In the next step, we analyze the quadruple integral 

(7.7) 



■1 



dz'dw'dzdw 



(2vri) 4 J ro J Tl 2 J t , q a e L{F 2 {z)+F 1 {z')) ( w - z )( w - z ')( w > _ z )(w> _ z >) ' 



which is the same integral as in (7.1) but now with deformed contours. 



Since the deformed contours are of steep descent /ascent and the functions 
Fj appear both the numerator and denominator, the integrand is bounded. 
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In fact, the main contribution, comes from small neighborhoods around the 
saddle points. As in Section 6.3, we introduce the local variables as in (6.6) 
around the saddle points. For example, near and O2 we introduce 

s t 

U, = n i + L l/2(_ jF v/ (ni)) l/2' Z = ^2+ L 1/2 (F // (5^)1/2' 

(7 ' 8) , S' , t' 

W ' ' = ^ 2+ LV2(_^/(n 2 ))l/2' Z, = 0l + L l/2 (jF // (fll)) l/2- 

The scaling by L 1 / 2 in the local variable, implies that the integrals tend to 
zero as L — > 00 (as was the case for the double integral for Kg in (6.27)). 
To this end, note that 

dwdzdw'dz' ±1 1 



(w -z)(w- z'){w' - z)(w' - z') LF{'(n)F"(n 2 ) (s ± it')(s> ± it) 

1 



fii - ^2 + 



L l/2(_ F »(Q 1 ))l/2 Ll/2(i^'(n 2 ))l/2 
1 

X ; 

2 - f2i + 



L 1 /2(F»(n 2 ))i/2 L 1 /2(^'(n 1 ))i/2 

The sign depends on the choice of the square roots, but is irrelevant in 
this discussion. Since by assumption the points f^i and Q2 are sufficiently 
far apart and the fact that the right-hand side is integrable (despite the 
singularities), we indeed see that the quadruple integral is of order C(L _1 ) 
as L — > 00. Moreover, the convergence is uniform on compact subsets of T>. 

We now come to the triple integrals presented in Figure 13 . Let us start 
with the top right case. In that situation the triple integral is given by 

rfi 2 e LF t (w) dzdz'dw 

(7.10) ' ' ' 



r Jr[ i2 Jn 2 e^(*0 (w -z)(w- z')(z - z<) 

The integral over the path from O2 to ^2 can be explicitly computed. As 
for the integrals over To and 2 , the main contribution comes again from 
small neighborhoods near the saddle points. By introducing local variables 
as in (7.10) and arguing as before, it is apparent that the integral vanishes 
as L — > 00 uniformly on compact subsets of T>. The same can be done for 
the situation in the lower right picture. 



For the triple integral over contours as in the top left picture of Figure 13 
we have 

rm e L(F 1 (w)+F 2 (z')) dz'dzdw 

(7.11) 



1,2 Jtia 



e L{ Fl {z')+F 2 (z)) [w - z )(w - z >)(z 
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In this triple integral we deform the path from rj 2 to 772 to an arc of circle 
centered at the origin which is path of steep descent for Re(i<2 — F%) for 



the case m\ < 771,2 (as mentioned earlier in the treatment of (7.5). Hence, 
it is exponentially small as L — > 00. Note that in the latter deformation 
we possibly pick up a residue again, which results in an additional double 
integral. However, since the contours are paths of steep descent /ascent the 
contribution of this term, possibly after an additional deformation, can be 
shown to be negligible. 

Concluding, from the seven multiple integrals obtained after deforming 
the contours, the leading term in the asymptotic expansion for R comes 
from the double integral (7.4). This integral equals the right-hand side of 



(4.17) and we proved the statement. 



Finally, we come to (4.18) 



this bound follows from (4.17) 



If (xj,m,j) are sufficiently far apart, then 
If on the other hand the points are close, 
then we need to deal with singularities in the integrals above. For example, 



the double integral that is over the contours in the left picture of Figure 12 
is no longer exponentially small. In fact, both double integrals (7.4) and 
(7.5) grow logarithmically with L as L — > 00, in case the points fL come 



closer and closer. The fact of the matter is that in all integrals obtained 
after deformation, the logarithmic behavior is however the worst that can 



happen. This can be checked by introducing local variables (7.10) in the 



triple and quadruple integrals. For example, in the quadruple integral (7.7) 



the integrand has a singular term (7.9). As long as il\ 7^ O2, then this is 
still integrable and the result is of order 0(log(f2i — Q2)) as Q% — CI2 — > 0. 



Hence for mi < 777.2 we obtain (4.18). 

The other situations 7714 = 7772 and m\ > 777-2 can be dealt with in a 
similar way. However, care should be taken in case (x\, m\) = (X2, 777,2) (and 
hence Oi = ^2)- In that case, deforming the contours leads to divergent 
integrals. The way around this is to perturb one of the saddle points (and 
the paths of steep descent /ascent) with a term of 0(L~ l ) and argue as 
above. Then we again obtain the logarithmic growth as L — > 00 and this 
shows that (4.18) holds for all points (xj,mj) in the bulk of LV. □ 
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